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YbRh 2 Si2: Quantum tricritical behavior in itinerant electron systems 
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We propose that proximity of the first-order transition manifested by the quantum tricrit- 
ical point (QTCP) explains non- Fermi-liquid properties of YbRli2Si2. Here, at the QTCP, a 
continuous phase transition changes into first order at zero temperature. The non-Fermi-liquid 
behaviors of YbRii2Si2 are puzzling in two aspects; diverging ferromagnetic susceptibility at 
the antiferromagnetic transition and unconventional power-law dependence in thermodynamic 
quantities. These puzzles are solved by an unconventional criticality derived from our spin 
fluctuation theory for the QTCP. 

KEYWORDS: quantum critical phenomena, quantum tricritical point, YbRh2Si2, non-Fermi-liquid behav- 
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Critical temperatures of the symmetry-breaking phase 
transitions can be lowered to zero at the quantum crit- 
ical point (QCP) by tuning quantum fluctuations such 
as by magnetic fields as shown in Fig. 1(a). Quantum 
critical phenomena in metals have attracted much in- 
terest from both theoretical and experimental points of 
view, because of not only its own right but also uncon- 
ventional superconductivity as well as non-Fermi-liquid 
behavior observed near the QCP. 1 ' 2 

The conventional spin fluctuation theory of the QCP 
by Moriya, Hertz and Millis 2-5 has succeeded in explain- 
ing a number of non-Fermi-liquid properties. However, 
this picture has been challenged by many recent ex- 
periments, 1 ' 2,6 where criticalities of thermodynamic and 
transport properties do not follow it. 

A typical heavy- fermion compound YbRIi2Si2 1 be- 
longs to such an unconventional category. At the mag- 
netic field H = 0, it exhibits an antiferromagnetic (AF) 
transition at the Neel temperature Tn = 0.07K. An AF 
QCP emerges at the critical magnetic field H c ~ 0.06T 
along the c axis. 7 ' 8 Near H c , Sommerfeld coefficient of 
specific heat 7 is logarithmically increased with lower- 
ing temperature T above 0.3K and even faster below 
it 9 in contrast to the conventional theory predicting 
convergence to a constant. Transport and optical data 
roughly show the resistivity linearly scaled with T and 
frequency. 7 Among all, a key aspect is an unusually en- 
hanced ferromagnetic susceptibility xo roughly scaled by 
Xo oc T-< and \o oc \H - ff c |~ c ' with C ~ C ~ 0-6 8 
contradicting the standard expectation of saturation to a 
constant. In accordance, the magnetization shows convex 
dependence on H. 8 NMR 10 and ESR 11 signals are also 
consistent. These non-Fermi-liquid properties are all con- 
tradicting the standard theory 2-5 for the AF QCP and 
are under extensive debates. 6 

A hint comes from the fact that the first-order transi- 
tion is observed for YbRli2Si2 under pressure. 12 Actually, 
the proximity of the first-order transition is common in 
many compounds with unconventional non-Fermi-liquid 
properties. Our idea is that the proximity of the first- 
order transition, namely, tricriticality solves the puzzle 
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Fig. 1. (color online), (a) Phase diagram of antiferromagnetic 
(AF) phase with critical line [solid (green) curve] ending at the 
QCP [(yellow) diamond] in T-H plane, where T (H) represents 
temperature (magnetic field), (b) Phase diagram with TCP [(yel- 
low) circle] separating the continuous [thin (green) curve] and 
first-order [thick (red) curve] transition lines, (c) Global phase 
diagram with tricritical line (TCL) separating the surfaces of 
continuous [above TCL (green)] and the first-order [below TCL 
(red)] surfaces. Here g represents parameter to control quan- 
tum fluctuations. In YbRh2Si2, g may correspond to pressure 
measured from the ambient one. The QTCP (circle) appears at 
(g,H,T) = (gt, Ht,0), namely the endpoint of TCL. Cross sec- 
tional view at constant g for g < gt[g > gt] corresponds to the 
phase diagram (a) [(b)]. 



because the tricriticality necessarily induces ferromag- 
netic tendency even at a clear AF transition. 

At T ^ 0, the tricritical point (TCP) where phase 
transitions change from continuous to first order as in 
Fig. 1 (b) has been studied in detail. 13 A characteristic 
feature of TCP is the diverging susceptibility not only at 
the ordering wavenumber Q but also at zero (xo)- 

If quantum fluctuations suppress the temperature of 
TCP to zero, quantum tricritical point (QTCP) appears 
[see Fig. 1(c)]. Then QTCP may alter the criticality of 
QCP as a proximity of the first-order transition. 

Recently TCP has been studied for the itinerant fer- 
romagnet 15 to understand the nature of the global phase 
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diagram of weak itinerant ferromagnets ZrZn2 and UGe2- 
Furthermore, ferromagnetic QTCP has been studied for 
itinerant helical ferromagnet MnSi 16 and nearly ferro- 
magnetic metal Sr3Ru207. 17 However, these previous 
studies on the ferromagnetic QTCP do not explain the 
unconventional coexistence of the ferromagnetic and AF 
fluctuations observed near the AF QCP in YbRh 2 Si2- 

In this letter, we propose that the proximity of the 
first-order transition opens a way to solve the puzzles 
in the AF quantum critical phenomena. The proximity 
of the first-order transition inherently generates diverg- 
ing ferromagnetic fluctuations concomitantly with the 
order-parameter (AF) fluctuations. The emergence of the 
concomitance is manifested by the quantum tricritical- 
ity, which generates an unexplored non-Fermi liquid. An 
unconventional scaling is derived by extending the self- 
consistent renormalization (SCR) theory 3 for spin fluc- 
tuations. Our result accounts for the otherwise puzzling 
properties of YbRli2Si2, even when we do not consider 
the possible valence transitions 18 or collapse of /-electron 
itinerancy as in the picture of the local quantum critical- 
ity- 19 

To understand the QTCP, we start from a standard 
Ginzburg-Landau- Wilson (GLW) expansion effective ac- 
tion for bozonic spin fields <p q at the wave number q: 3 ~ 5 



E 

q,q',q" 



u(q, q\q")(f q ■ <p- q ') 



between M and Aft as 

AT = a + aiM f2 + a 2 M^ + • • • , (5) 

where the expansion coefficients clq ~ a 2 are determined 
by substituting eq. (5) into the saddle point condition: 

f (T, H)a + 4u (T, /f )a 3 + Gva% - H = 0, (6) 

Ua u Q (T,H) + a 1 R(T,H) = 0, (7) 

where R(T, H) = f (T,H) + l2uo(T,H)a 2 + 30va%. By 
using eq. (5), we obtain the free energy as 

F =\r Q {T, H)M^ + u> Q (T, H)M^ + 0(M+ 6 ), (8) 

where u' Q (T,H) = u Q (T,H)(l + 6a Oi). In eq. (8), 
continuous phase transitions occur at tq — when 
u'q(T,H) > 0, while the first-order phase transitions oc- 
cur when u' Q (T,H) < 0. 13 ' 14 Therefore, the QTCP ap- 
pears when the conditions fg(0, H t ) — and uq(0, H t ) = 
are both satisfied, 20 where H t is the critical field at the 
QTCP. 

We now discuss the susceptibilities xq at the AF vec- 
tor Q and xo at q = in the disordered phase (Aft = 0, 
M = ao) by using eq. (6) and the free energy (8). From 
eq- (8), Xq 1 is given as 

= f Q (T,H). (9) 



Xq 



Mt=0 



qi ' ( P-q2 



By differentiating eq. (6) with respect to the magnetic 
field H, we obts : ~ 



am Xo 
l 



as 



9l~95 



Xo 



(1) 



where H is external magnetic field; u(q, q', q") and v are 
constants, while r q depends on the magnetic field H. 
From eq. (1), the free energy F is obtained from 



V OH ) 1 - a dr /dH - Aa^duo/dH 
<xuq(T,H). (10) 

Here, we used eq. (7), which gives R(T, H) oc uq(T, H). 
Now, the fluctuation-dissipation (FD) theorem 21 



exp(-F/T) = fj]Vip q exp(-S[ip q ]/T). (2) E^l") = - / du{- + n(w)J ^ Im X (g,w), (11) 

J „ q^0,Q J0 9#0,Q 



Since the QTCP is expressed by fluctuations at both the 
AF Bragg wavenumber Q and 0, we approximate the free 
energy as a function of the order parameter Aft = (<pq) 
and the uniform magnetization M = (fo)- 



F =^f Q M^ 2 + u Q M^ + wM t6 



1 



-f M 2 + w AP + vM 6 - HM, 



(3) 



where ?q, uq, fo, uo and K, are defined as fg(T, H) = 
r Q (H) + 12u Q (IC + M 2 ) + 90v(!C + M 2 ) 2 , u Q (T,H) = 
u Q + 15v(IC + A/ 2 ), f (T, H) = r (H) + 12u JC + 90v)C 2 , 
uq(T, H) = Mo + 15vJC, and 



■■ E (i^i 2 )- 

9#0,Q 



(4) 



Effects of spin fluctuations are included in /C following 
the SCR theory. We approximate u(q,q,Q) [u(q, q, 0)] 
and the equivalent coefficients as q- independent values; 
u(q, q, Q) ~ u Q [u(q, q, 0) ~ u ] for all q [for q ^ Q] . 

We eliminate M in eq. (3) by using the saddle point 
condition for AT, 8F /dM = leading to the relation 



and n(u>) = l/(e"/ T — 1), combined with cqs. (4), (6), 
(9), and (10) constitute the self-consistent equations 3 to 
determine K. and x in the scheme of our extended SCR 
theory. Using this SCR theory, we now clarify how the 
susceptibilities and the magnetization measured from the 
QTCP (x Q \xo\Sao — ao — aot with aot being the value 
at the QTCP) are scaled with 5H = H — H t and T near 
the QTCP. The results will be shown in cqs. (15)-(17). 

In the SCR theory, non-trivial temperature depen- 
dence of physical properties comes from the spin fluc- 
tuation term JC. Therefore, we first clarify the scaling 
of K, by using the FD theorem combined with expan- 
sions of Xo+q(^) and XQ+q{ u ) i n terms of the wavenum- 
ber q and the frequency uo near the QTCP. The order- 
ing susceptibility XQ+q{u) is assumed to follow the con- 
ventional Ornstein-Zernike form (xq+ 9 (^) _1 — Xq 1 + 
AQq 2 — iCqlo) as in the SCR formalism, while the uni- 
form part xo+q{u) turns out not to follow. This is be- 

— 1 — 1/2 

cause the scaling relation Xo ^ Xq holds near TCP 
within the GLW theory. 13 As we will see, the self- 
consistency among eqs. (4), (6), (9), (10), and (11) re- 
quires that this relation still holds for q ^ 0. Therefore, 
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we obtain the relation as Xo+<?(0) 1 oc x<3+<?(0) 1 I 2 oc 
(Xq 1 + A Q q 2 ) 1/2 oc (xq 2 + A q 2 ) 1/2 . From the conserva- 



oc 



T 3/2 



tion law, ui dependence of xo+«j( w ) 1 should be given 

1 1 as xo+gi^y 1 



as Xo+ q (u) 1 - Xo+g(0) 1 - iC u;/q. Finally, we ob- 



tain u> and q expansions of xo+gO^O 
(Xo 2 + A q 2 ) 1/2 ~^/q. ' 

By substituting the above forms for xo+q(v) 1 
[XQ+qiw)^ 1 } into the FD theorem (11), the contributions 
from the spin fluctuations near zero [ordering] wavenum- 
ber is given as E r o(^) - K oo - K oiX 2 + K 0T T 2 , 
E,~q<^> - K Q0 ~ KqiXq 1 + K QT T 3 / 2 } where K 00 , 
if oi, and Kqt [Kqq, i^Qii and Kqt] are constants. From 
these relations, in three dimensions, we obtain the scaling 
of SIC measured from the QTCP as 

SK ~ -K 01 Xo 2 - KqiXq 1 + KotT 2 + K QT T 3 / 2 . (12) 

The singularity of magnetization ao is obtained by 
solving eq. (6). Near the QTCP, eq. (6) can be approxi- 
mated as ASa 2 , + BSa + C — 0, with A — 12a t(5val t + 
tto), B = Sfo + 12a,Q t 5uQ, and C = aotSfo + 4ap t (5uo — SH, 
where Sfg = fo(T, H) — fo(0, H t ), and Suq — uq(T, H) — 
uo(0,H t ). Since both B and C vanish at the QTCP, we 
obtain the asymptotic behavior of Sao as 

(a SH + ai SIC) 1/2 , 



5a Q 



(13) 



where ao,ai are constants. 

By defining 5f Q (T, H) = r Q (T, H) - r Q (0, H t ), we get 

Xq 1 = 6f Q (T, H) =Sr Q (H) + 90v(SIC + <5a ) 2 , (14) 

since both fq(0, H t ) and uq(0, H t ) are zero at the QTCP 
and terms linear in SK, and Sda vanish. Here SrQ and 
5ao are defined as Stq = tq{H) — rg(ift) ~ tqh5H, 
8ao = Oq — Oqi = 5a {5ao + 2a t)- 

From eqs. (12) and (13), SK. is higher order of Sao near 
the QTCP. Then, from eqs. (10) and (14), the most domi- 
nant terms of Xq 1 and x^ 1 are given as Xo" 1 00 ^ a o, Xq 1 — 
tqhSH + 360val t Sa^, which together with eqs. (12) and 
(13) leads to SH and T dependence as 

Xq 1 ~f3 QQ SH + P Q1 T 3 / 2 , (15) 
Xo^^oo^ + ^oiT 3 / 2 ) 1 / 2 , (16) 
Sa - (a SH + a[T 3/2 ) 1/2 . (17) 

Singularities of the uniform susceptibility xo, the order- 
ing susceptibility XQi and the magnetization Sao near 
the QTCP are summarized in Fig. 2. We note that the 
singularity of Xq 1 given by Green et al 17 as T 8 / 3 is not 
correct, since they neglect the T 2 dependence of the bare 
second-order coefficient r q . 5 

We now examine whether the criticality of the QTCP 
is consistent with the experimental results of YbPdi2Si2. 
We first emphasize that the puzzling critical exponents 
in YbRh 2 Si2 described above for T and H dependences 
of Xo and M are well consistent with the quantum tricrit- 
icality derived from eqs. (16) and (17), while any other 
theories do not reproduce these exponents. 

To compare with the experimental results more quan- 
titatively, we now solve the self-consistent equations 
numerically, and obtain the uniform susceptibility xo 
and the magnetization Sao near the QTCP as follows: 




5M oc T 3 / 4 x -l oc SH 

Xo" 1 SH 1 / 2 
5M oc 5H l l 2 



Fig. 2. (color online). Schematic phase diagram around the 
QTCP under the magnetic fields H. 8M denotes the magne- 
tization measured from the critical value at the QTCP. 



We first approximate the magnetic field dependence of 
Sr (H) = r (H) - r (H t ) as Sr Q (H) ~ r 0H SH. The con- 
tributions from the spin fluctuations, namely /C, can be 
calculated by setting the four SCR parameters T a, Too, 
Tqa, and Tq 22 Furthermore, once the parameters v, 
r QH, toh, Ht, and dot are fixed, the other parameters 
(ro, tq, ito, and uq) are determined from the condi- 
tions r Q (0,H t ) = 0, u Q {0,H t ) = and eqs. (6), (7). 
In this letter, to calculate physical properties, we em- 
ploy a reasonable set of parameters given in ref. 23 as 
follows: We estimate H t and dot directly from exper- 
iments and also choose conventional SCR parameters 
(To a , Too j Tqa , Tqq ) within the order of 10-100K. This 
range of SCR parameters is typical in heavy fermion 
compounds. 3 In contrast to these, for the other non- 
primary parameters {tqh, ran, and v), we do not find 
any constraint from physical requirement. Therefore, we 
have freely tuned these parameters to reproduce the ex- 
perimental results quantitatively. However, the critical 
exponents do not change even if we have chosen these 
parameters arbitrarily. Microscopic derivation of these 
phenomelogical parameters is left for future studies. 

In Fig. 3 (a), the numerical result of the temperature 
dependence of Xo" 1 j us t on t ne QTCP is compared with 
the experimental Xo" 1 reported in ref. 8. Although we 
obtain the critical exponent £ = 0.75 (xo 1 k T^) for 
T — > 0, the numerical result shows that Xo 1 lS roughly 
scaled by T 6 at higher temperatures (T >1.0K). We 
emphasize that the puzzling convex behavior of Xo 1 
for YbRli2(Sio.95Geo.o5)2 near the QCP can not be ac- 
counted for by the conventional quantum criticality, be- 
cause the critical exponent Q is always larger than one for 
the conventional quantum critical point. 3-5 The nonzero 
offset of experimental Xo" 1 at T = indicates that the 
QCP in YbRh 2 (Sio.95Geo.o5)2 exists slightly away from 
the QTCP. It is an intriguing experimental challenge to 
determine the precise location of the QTCP by tuning 
the pressure and the magnetic field. 

In Fig. 3 (b), the numerical result of magnetization 
curve at T = is compared with the experimental re- 
sult reported in ref. 8. By using the same parameters, 
it is noteworthy that not only the uniform susceptibil- 
ity but also the magnetization is consistent with the ex- 
perimental result quantitatively in the relevant parame- 
ter region. Because the experimentally observed QCP is 
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Fig. 3. (color online). (a) Experimental Xo 1 f° r 
YbRh 2 (Sio.95Geo.05)2 at H = 0.03(T) reported in ref. 8 
compared with the present SCR theory. Green line (red circle) 
represents the experimental (theoretical) \o Theoretical Xo" 1 
is calculated just on the QTCP (H = Ht). (b) Experimental 
magnetization curve for YbRh2(Sio.95Geo.os)2 at T = 0.09(K) 
reported in ref. 8 compared with the present theory. Green 
line (red circle) represents the experimental (theoretical) 
magnetization curve. 5M (SH) represents the magnetization 
(magnetic field) measured from the critical value. We estimate 
the experimental critical magnetic field H c (magnetization M c ) 
as 0.027(T) (0.004 (/ib))- The inset shows the numerical result 
of Sommerfeld coefficient of the specific heat 7 just on the 
QTCP. 



slightly away from the QTCP, small deviations are seen 
at low temperatures and low magnetic fields (T < f .OK 
and 5H < 0.5T). We note that this singularity of the 
magnetization [SM oc 5H 1 / 2 as seen from eq. (17) and 
Fig. 3 (b)] is qualitatively different from that of the con- 
ventional metamagnetic transitions, ft has been proposed 
that they belong to the Ising universality class. 24 The 
critical exponent 5 (SM oc SH 1 / 6 ) of the Ising universal- 
ity is always larger than three at any dimensions. The 
present critical exponent 5 = 2 makes a sharp contrast 
to such conventional critical exponents 5 > 3. 

We now discuss the singularity of the Sommerfeld co- 
efficient of the specific heat 7. In the inset of Fig. 3 (b), 
the numerical result of 7 just on the QTCP is shown. 
At high temperatures (T > 1.0K), both the singular- 
ity (7 oc — log T) and the amplitude are consistent with 
those of the experimental result. 9 However, at low tem- 
peratures (T < 1.0K), within this SCR theory, 7 near 
the AF QTCP approaches that of the conventional AF 
QCP (7 oc const. — T 1 / 2 ), while experimentally, power- 
law-like behavior is observed for T < 0.37^. 9 This dis- 
crepancy may be solved by considering either the fact 
that the Neel temperature is actually nonzero or effects 
of valence fluctuations, 18 while the criticality of magnetic 
properties analyzed here should remain unchanged. 

Finally, we discuss a different scenario for the QCP 
in YbRli2Si2 proposed by Coleman et al. 19 They claim 
that a breakdown of composite heavy fermion (namely, 
all / electrons decouple from the Fermi surface) occurs 
at the QCP and no heavy electron exists in the ordered 
phase any more. This scenario is inconsistent with a large 
Sommerfeld coefficient of the specific heat 7 even in the 
ordered phase. 9 While a large change of Hall constant 
in YbRh2Si2 25 was suggested to support their scenario, 
Norman 26 has pointed out that small changes of the / 



electron occupation are sufficient to reproduce the ex- 
perimental result by calculating the band structures of 
YbRh2Si2. Steep change in the Hall coefficient is then 
naturally understood under the proximity of the first- 
order transition. 

In summary, a non-Fcrmi liquid different from that 
obtained from the conventional QCP is shown to emerge 
when the proximity of the first-order transition is in- 
volved through the QTCP. The unconventional criticality 
thus obtained by the extension of the SCR theory solves 
the puzzles in the experimental results of YbRli2Si2. It is 
intriguing to examine whether this proximity also plays 
roles in other unconventional non-Fermi liquids. 
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